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ABSTRACT 

We compare optical and hard X-ray identifications of active galactic nuclei (AGNs) using a uni- 
formly selected (above a flux limit of /2_s kov = 3.5 x ICT 15 erg cm~ 2 s _1 ) and highly optically 
spectroscopically complete (> 80% for / 2 _s kcV > 10 -14 erg cm -2 s _1 and > 60% below) 2 — 8 keV 
sample observed in three Chandra fields (CLANS, CLASXS, and the CDF-N). We find that empirical 
emission-line ratio diagnostic diagrams misidentify 20 — 50% of the X-ray selected AGNs that can be 
put on these diagrams as star formers, depending on which division is used. We confirm that there 
is a large (2 orders in magnitude) dispersion in the log ratio of the [OIII]A5007 (hereafter, [OIII]) to 
hard X-ray luminosities for the non-broad line AGNs, even after applying reddening corrections to the 
[OIII] luminosities. We find that the dispersion is similar for the broad-line AGNs, where there is not 
expected to be much X-ray absorption from an obscuring torus around the AGN nor much obscuration 
from the galaxy along the line-of-sight if the AGN is aligned with the galaxy. We postulate that the 
X-ray selected AGNs that are misidentified by the diagnostic diagrams have low [OIII] luminosities 
due to the complexity of the structure of the narrow-line region, which causes many ionizing photons 
from the AGN not to be absorbed. This would mean that the [OIII] luminosity can only be used to 
predict the X-ray luminosity to within a factor of ~ 3 (one sigma). Despite selection effects, we show 
that the shapes and normalizations of the [OIII] and transformed hard X-ray luminosity functions 
show reasonable agreement, suggesting that the [OIII] samples are not finding substantially more 
AGNs at low redshifts than hard X-ray samples. 

Subject headings: cosmology: observations — galaxies: active — galaxies: nuclei — galaxies: Seyfert 
— galaxies: distances and redshifts — X-rays: galaxies 



1. INTRODUCTION 

Determining which sources in a large sample of galax- 
ies host active galactic nuclei (AGNs) at their centers 
is a challenging but essential step for galaxy evolution 
studies. The most common method used optically to 
separate star-forming galaxies from AGNs is a Baldwin 
et al. (1981)-type empirical diagnostic diagram (here- 
after, BPT diagram; see also Osterbrock & Pogge 1985 
and Veilleux & Osterbrock 1987) of the narrow emission- 
line ratios [OIII] A5007/H/3 versus [NII]A6584/Hx In this 
diagram galaxies occupy two well-defined wings: the left 
wing consists of star-forming galaxies, while the right 
wing is attributed to galaxies with an active nucleus (e.g., 
Kauffmann et al. 2003; Stasihska et al. 2006; Kewley 
et al. 2006). The basic idea underlying these diagrams is 
that the emission lines in star-forming galaxies are pow- 
ered by massive stars, so there is a well-defined upper 
limit on the intensities of the collisionally excited lines 
relative to the recombination lines (such as Ha or H/3). 
In contrast, AGNs are powered by a source of far more 
energetic photons, making the collisionally excited lines 
more intense relative to the recombination lines. 

1 Some of the data presented herein were obtained at the W. 
M. Keck Observatory, which is operated as a scientific partner- 
ship among the California Institute of Technology, the University 
of California, and the National Aeronautics and Space Administra- 
tion. The observatory was made possible by the generous financial 
support of the W. M. Keck Foundation. 

Department of Astronomy, University of Wisconsin-Madison, 
475 N. Charter Street, Madison, WI 53706 

3 Department of Physics and Astronomy, University of Hawaii, 
2505 Correa Road, Honolulu, HI 96822 

4 Institute for Astronomy, University of Hawaii, 2680 Woodlawn 
Drive, Honolulu, HI 96822 



Kauffmann et al. (2003) and Heckman et al. (2004) 
proposed using the luminosity of the [OIII]A5007 (here- 
after, [OIII]) line as a tracer of AGN activity. As they 
pointed out, the advantages of the [OIII] line are that (1) 
it is typically strong and easy to detect; and (2) although 
it can be excited by both massive stars and AGNs, it has 
been observed to be relatively weak in metal-rich, star- 
forming galaxies. However, the disadvantage is that one 
needs to assume that [OIII] is an unbiased, orientation- 
independent indicator of the ionizing flux from the AGN. 

In the standard 'unified' model for AGNs (e.g., An- 
tonucci 1993), if the observer's view of a galaxy's cen- 
tral supermassive black hole and its associated contin- 
uum and broad emission-line region are unobscured (ob- 
scured) by the presence of a dusty torus, then it is a 
type 1 (type 2) AGN. However, even in the type 2 AGNs 
the covering factor is not one, so radiation is able to es- 
cape through the opening angle and photoionize the gas 
located several hundred parsecs or more from the cen- 
tral engine. This gives rise to the narrow emission-line 
region with its strong narrow permitted and forbidden 
emission lines. Since the narrow emission-line region lies 
outside of the dusty torus, the emission lines should not 
suffer from obscuration by that high column density ma- 
terial (though they may be affected by dust within the 
host galaxy). Although there is some evidence that the 
[OIII] luminosity is a good measure of AGN activity (e.g., 
Mulchaey et al. 1994), and, indeed, many researchers 
have assumed it to be so in their analyses (e.g., Alonso- 
Herrero et al. 1997; Hao et al. 2005; Netzer et al. 2006; 
Bongiorno et al. 2010), recent work has called this con- 
clusion into question (e.g., Cocchia et al. 2007; Melendez 
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et al. 2008; Diamond-Stanic et al. 2009; La Massa et al. 
2009). 

With the advent of the Chandra and XMM-Newton X- 
ray observatories, AGN activity can now also be traced 
through the 2 — 10 keV luminosity, although this measure- 
ment may suffer from some absorption (e.g., Diamond- 
Stanic et al. 2009; La Massa et al. 2009). Also, X-ray se- 
lection misses the most heavily absorbed, Compton-thick 
sources with Nh > 10 24 cm~ 2 (see Comastri 2004 for a 
review). An interesting question is how well the X-ray 
and optical selections compare. With our large, uniform, 
and highly optically spectroscopically complete OPTX 
sample of X-ray selected AGNs, we are in a good posi- 
tion to do this comparison. The OPTX sample consists 
of one deep pencil-beam survey (Chandra Deep Field- 
North or CDF-N) and two moderately deep wide-field 
surveys ( Chandra Large Area Synoptic X-ray Survey or 
CLASXS and Chandra Lockman Area North Survey or 
CLANS). 

In this paper we explore how well the BPT dia- 
gram does at classifying as AGNs the sources in the 
OPTX X-ray selected AGN sample and how that clas- 
sification relates to the very broad observed range in 
£2-10 kev/£[om] ( e -g-j Cocchia et al. 2007; Melendez 
et al. 2008; La Massa et al. 2009). Following Heck- 
man et al. (2005), Bongiorno et al. (2010), and Geor- 
gantopoulos & Akylas (2010), we use our measurements 
of the mean luminosity ratios to transform the fits to 
the OPTX+ SWIFT BAT hard X-ray luminosity func- 
tion (LF) from Yencho et al. (2009) and the local RXTE 
hard X-ray LF from Sazonov & Revnivtsev (2004) into 
[OIII] LFs, which we then compare with the SDSS [OIII] 
LFs from Hao et al. (2005) and Reyes et al. (2008) to in- 
vestigate the efficiency of optical versus X-ray selection. 

The structure of the paper is as follows. In Sec- 
tion 2 we briefly describe our 2 — 8 keV Chandra selected 
OPTX sample. In Section 3 we determine optical nar- 
row emission-line luminosities and investigate the effects 
of line reddening. In Section 4 we test how many X- 
ray selected AGNs are misidentificd as star formers by 
the BPT diagram. In Section 5 we plot the ratio of the 
[OIII] luminosities to the 2 — 8 keV luminosities versus 
the 2 — 8 keV luminosities for our non-BLAGNs and look 
for any luminosity dependence, as suggested by Netzer 
et al. (2006). In Section 6 we discuss the probable impact 
of the observed dispersion in this luminosity ratio on the 
BPT classification of X-ray AGNs. We then align various 
hard X-ray and [OIII] LFs for comparison. In Section 7 
we summarize our results and their implications. 

All magnitudes are in the AB magnitude system, and 
we assume VLm = 0.3, = 0.7, and H = 70 km s _1 
Mpc' 1 . 

2. SAMPLE 

In Trouille et al. (2008) we presented the X-ray cata- 
log for the CLANS field. We also presented the optical 
and infrared photometry and the spectroscopic and pho- 
tometric redshifts for the X-ray sources in all three of 
the OPTX fields (CLANS, CLASXS, and CDF-N). In 
Trouille et al. (2009) we updated the X-ray catalog for 
the CLANS field, including 28 sources that had inad- 
vertently been dropped in the original catalog, and we 
provided additional spectroscopic redshifts for this field. 

In Trouille et al. (2009) we used a uniform and highly 



significant subsample of the OPTX sources to compare 
the optical spectral types with the X-ray spectral prop- 
erties. Here we use the same selection of sources in the 
2 — 8 keV band with significance greater than 3 a and 
fluxes greater than the flux limits (for a S/N — 3) at the 
pointing center of the 70 ks CLANS pointings. There 
are 746 X-ray sources in this 2 — 8 keV sample (410, 
251, and 85 sources from the CLANS, CLASXS, and 
CDF-N fields, respectively). Our optical spectroscopic 
completeness with the DEep Imaging Multi-Object Spec- 
trograph (DEIMOS; Faber et al. 2003) on the Keck II 
10 m telescope is very high (see Table 3 of Trouille 
et al. 2009), particularly at the bright flux end. At 
/2-8 kov > 1-4 x 10~ 14 erg cm~ 2 s _1 (the break flux 
in the 2 — 8 keV number counts; see Trouille et al. 2008) 
we have spectroscopic redshifts for 178 of the 217 sources 
(i.e., we are 82% complete to this flux). 

In Trouille et al. (2009) we followed Szokoly et al. 
(2004) and Barger et al. (2005) and classified the spec- 
troscopically identified sources into four optical spectral 
types: absorbers (ABS; no strong emission lines); (2) 
star formers (SF; strong Balmer lines and no broad or 
high- ionization lines); (3) high-excitation sources (HEX; 
[NeV], CIV, narrow Mgll lines, or strong [OIII]); and (4) 
broad-line AGNs (BLAGNs; optical lines having FWHM 
line widths > 2000 km s 1 ). Although the HEX spectral 
type largely overlaps the classical Seyfert 2 spectral type, 
describing the sources as HEX sources helps to avoid con- 
fusion with the classical definitions. 

In Trouille et al. (2008, 2009) we used the X-ray fluxes 
and spectroscopic plus photometric redshifts to calculate 
the rest-frame 2 — 8 keV luminosities, Lx- At z < 3 
(which is all we will be considering in this paper), we 
calculated the luminosities from the observed-frame 2 — 
8 keV fluxes, assuming an intrinsic T — 1.8. That is, 
Lx = / x 47rd| x K— correction, where for z < 3 

K - corr = (1 + z)" ' 2 and / = / 2 _ 8 kcV . (1) 

Note that using the individual photon indices (rather 
than the universal power-law index of V = 1.8 adopted 
here) to calculate the K~ corrections would result in only 
a small difference in the rest-frame luminosities (Barger 
et al. 2002). We have not corrected the X-ray luminosi- 
ties for absorption. We will examine the impact of this 
approach in Section 6. 

Any source more luminous than Lx — 10 42 erg s _1 is 
very likely to be an AGN on energetic grounds (Zczas 
et al. 1998; Moran et al. 1999). Of the 746 sources in 
our 2 — 8 keV sample, 619 have Lx > 10 42 erg s _1 . If 
we adopt the Ranalli et al. (2003) relation between star- 
formation rate (SFR) and X-ray luminosity 

£2-10 kcV = SFR/M Q x 5 x 10 39 erg s" 1 , (2) 

we find that a SFR of at least 200 M yr _1 is required to 
produce an X-ray luminosity greater than 10 42 erg s _1 . 
Given that the space density of ultraluminous infrared 
galaxies (ULIRGs, which require a SFR > 200 M Q ) is 
approximately 10~ 5 Mpc~ 3 at z ~ 0.9 (Magnelli et al. 
2009), near the peak in our redshift distribution, and 
our OPTX fields cover <~ 1.2 deg 2 , we expect ~ 30 
ULIRGs total in our sample. Because the space den- 
sity of ULIRGs increases rapidly with redshift, if instead 
we restrict our analysis to z < 0.5 (see Section 4), we ex- 
pect fewer than three ULIRGs contaminating our sample. 
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Since the contamination rate is low 



pie), we assume that the 619 Lx > 10 42 erg s 1 



5% of our sam- 
sources 

are AGNs (although see Section 4 for our classifications 
of these sources as star formation or AGN dominated 
using emission-line ratio diagnostics). In the remainder 
of this work, our '2 — 8 keV sample' refers to these 619 
sources. 

3. HIGH-IONIZATION NARROW EMISSION LINES 
3.1. [OIH1A5007 

[OIII] appears in our optical spectra for sources with 
z < 0.95. Since DEIMOS degrades at the red end, we 
limit our [OIII] analysis to the 157 sources in our 2—8 keV 
sample with spectroscopic redshifts less than 0.85. We 
use the method described in Kakazu et al. (2007) to de- 
termine the [OIII] flux. We first determine which band- 
pass covers the redshifted emission line. We then inte- 
grate the spectrum convolved with that bandpass filter 
response and set the result equal to the broadband flux 
to determine the normalization factor. We determine 
the continuum using a sliding 250 A median. We then 
multiply the integrated continuum-subtracted emission 
line by the normalization factor to obtain the flux in the 
line. This procedure only works for sources with secure 
continuum magnitudes where the sky subtraction can be 
well determined in the spectra. We then convert the 
emission-line flux into an emission-line luminosity using 
£[Oiii] = /[oiii] x 47T(i|. 

In Figure 1 we show rest-frame 2 — 8 keV luminosity 
versus redshift for the sources in our 2 — 8 keV sample for 
which we were able to determine -Zj[oiii] (black symbols; 
122 of the 157 sources). We denote non-BLAGNs by tri- 
angles and BLAGNs by squares. Of the 35 sources that 
are lacking Lpm] determinations, 22 are classified as 
absorbers (magenta circles). By definition, absorbers ex- 
hibit little to no strength in emission lines. However, they 
do exhibit H/3 absorption (and Ha as well for those with 
a low enough redshift). The presence of these B aimer 
absorption lines suggests that they are not BL Lac ob- 
jects (see Laurent-Muehleisen et al. 1998 for a detailed 
description of BL Lac properties) . We note that if these 
are in fact low luminosity BL Lac objects, their low non- 
thermal optical luminosity could allow for the presence 
of stellar absorption lines. We use an equivalent width 
for the [OIII] line that is la above the noise to determine 
the upper limit on Lroiii] for these sources. 

The remaining 13 sources lack ironi] determinations 
as a result of observational complications. Six sources 
(red symbols) are fainter than the limiting magnitude in 
the filter needed to determine their Z/roiii] values. Four 
sources (green symbols) have sky lines that interfere with 
the ifoni] determinations, and three sources (blue sym- 
bols) have an intervening absorption line located at the 
position of [OIII]. We remove these 13 sources from our 
subsequent analysis. 

3.2. Emission- Line Reddening 

The average reddening towards the narrow-line re- 
gion in AGNs is still a matter of debate. Studies using 
HST spectroscopy combined with the IUE-b&sed anal- 
ysis of Ferland & Osterbrock (1986) suggest that the 
narrow emission-line reddening for non-BLAGNs corre- 
sponds to E(B — V) = 0.2 — 0.4 mag (e.g., Ferruit et al. 
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Fig. 1. — Log(Ljf ) vs. redshift for sources in our 2 — 8 keV sample 
(black — sources with £[om] determinations; magenta — absorbers; 
red — sources lacking magnitudes in the necessary filter; green — 
sources with sky line interference; blue — sources with intervening 
absorption interference). Non-BLAGNs are denoted by triangles, 
BLAGNs by squares, and absorbers by circles, as noted in the 
legend. The black curve corresponds to the flux limit for our 2 — 
8 keV sample. 
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Fig. 2. — (a) Ratio of Ha flux to H/3 flux vs. rest-frame 2 — 
8 keV luminosity for sources in our 2 — 8 keV sample with reliable 
/un and /h/3 measurements (dotted line — the canonical intrinsic 
Balmcr decrement (Jh //h^)o = 3.1; dashed line — the Zakamska 
et al. (2003) reddening value for their composite SDSS type 2 AGN 
spectrum, fua/fnp = 4.0; solid line — the Bassani et al. (1999) 
median reddening value for their sample of type 2 AGNs, /hc«//h^ 
= 6.4. (b) Same as (a) but versus redshift. 

1999). Zakamska et al. (2003) used the Balmer decre- 
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ment, fn a /fnp, measured from a composite spectrum of 
their SDSS non-BLAGNs as a proxy for the volume aver- 
aged narrow emission-line reddening and found a similar 
result of E(B — V) = 0.27 mag (though see Netzer et al. 
2006 for a discussion of possible selection effects). Dahari 
& De Robertis (1988) and Bassani et al. (1999) found sig- 
nificantly larger median reddening in their non-BLAGNs 
of E(B - V) = 0.48 mag and E(B - V) = 0.65 mag, 
respectively. The Dahari & De Robertis (1988) sam- 
ple represented the largest available uniform database 
of optical spectra at that time, including most of the 
known Markarian Seyferts and a majority of the classi- 
cal Seyferts listed in Weedman (1977, 1978). The me- 
dian reddening quoted above is for the 20 Seyfert Is in 
their sample with z < 0.03 and axial ratios b/a > 0.5. 
The Bassani et al. (1999) sample was composed of the 57 
Seyfert 2s with reliable hard X-ray spectra in the litera- 
ture at that time. They obtained the Balmer decrements 
from a variety of sources. The majority of their sources 
have z < 0.06, with a few at z < 0.2. 

In Figure 2 we show /h q //h/3 versus (a) rest-frame 
2 — 8 kcV luminosity and (b) rcdshift for the 16 sources 
in our 2 — 8 keV sample for which we were able to measure 
both the Ha and H/3 fluxes. For reference, there are 45 
non-BLAGNs and 19 BLAGNs with < z < 0.5 (the 
rcdshift at which Ha leaves our spectral window) in our 
2 — 8 keV sample. The separation between BLAGNs 
(gray squares) and non-BLAGNs (black triangles) with 
luminosity is a result of BLAGNs dominating the number 
densities at higher X-ray luminosities (Steffen et al. 2003; 
Ueda et al. 2003; Barger et al. 2005; La Franca et al. 2005; 
Simpson 2005; Akylas et al. 2006; Beckmann et al. 2006; 
Sazonov et al. 2007; Hasinger 2008; Silverman et al. 2008; 
Winter et al. 2009; Yencho et al. 2009). This result had 
previously been noted in the radio by Lawrence & Elvis 
(1982). 

To determine the Ha and H/3 fluxes for the non- 
BLAGNs we use the same procedure as described in Sec- 
tion 3.1 for determining the [OIII] emission-line fluxes. 
However, to properly account for stellar Balmer absorp- 
tion, we first fit a stellar population model to the contin- 
uum using the Tremonti et al. special-purpose template 
fitting code (private communication), described in detail 
in Tremonti et al. (2004) . In brief, we use a grid of stellar 
template models covering a range of age and metallicity 
generated by the Bruzual & Chariot (2003) population 
synthesis code. For each galaxy we transform the tem- 
plates to the appropriate rcdshift and velocity dispersion 
to match the data. We then subtract the best-fit model 
from the continuum. We only do the subtraction in our 
determination of the H/3 emission-line flux, because the 
Ha line is typically stronger with a relatively weaker un- 
derlying stellar absorption. We find an average equiva- 
lent width of <~ 1 A for the H/3 stellar absorption. This 
is in agreement with the fixed offset of 1 A used in Cowie 
& Barger (2008), which they determined from their av- 
eraged absorption-line galaxy spectrum. 

To determine the Ha and H/3 fluxes for the BLAGNs, 
we fit both a broad and a narrow component to the lines, 
subtract the broad component, and then follow the same 
procedure as described above. We are able to deter- 
mine reliably the Ha and H/3 fluxes for 26% (22%) of 
the BLAGNs (non-BLAGNs) in our 2 - 8 keV sample 
with < z < 0.5. 



We do not see any statistically significant difference 
in the reddening between our BLAGNs and our non- 
BLAGNs, nor do we see a dependence on redshift or 
X-ray luminosity. We note the small numbers in our sam- 
ple and the three non-BLAGNs with significantly higher 
Ma/flip values. The average /hq//h/3 for the com- 
bined samples is 5.2 ±3.1 (corresponding to E(B — V) = 
0.47±0.4), which falls between the Zakamska et al. (2003, 
dashed line) and Bassani et al. (1999, solid line) results. 
The dotted line shows the canonical intrinsic Balmer 
decrement value, (/h q //h/3)o = 3.1 (Ferland & Netzer 
1983; Osterbrock & Martel 1993). 

Following Bassani et al. (1999), wc correct our narrow 
emission line fluxes for optical reddening with the equa- 
tion 

/[OIII], corr = /[OIII] X [{ftta / fttp) / {ftta / fuf)) oY'^ ■ (3) 

We use our calculated /ho/ '/h/3 values for the 16 sources 
in which both lines are measured and the average 
fHu/.fnp= 5.2 for the remaining sources. This results 
in an average ratio between our extinction-corrected 
and non-extinction-corrected luminosities of 4.6. How- 
ever, we caution that there remain many uncertainties 
when applying standard extinction corrections to AGN 
emission-line luminosities. For example, Hao et al. (2005) 
argue that AGNs may have higher intrinsic fua/fap ra- 
tios than normal galaxies because of higher densities and 
radiative transfer effects. Also, the relative distribution 
of dust and emission-line gas does not necessarily follow 
the simple model in which dust forms a homogeneous 
screen in front of the line-emitting regions. 

A Kolmogorov-Smirnov (K-S) test reveals no signif- 
icant difference in the distributions of £[om] fo r our 
z < 0.5 BLAGNs and non-BLAGNs {p = 0.5). We 
do not include the upper limit values for our absorbers 
in this comparison. This is in contrast with Diamond- 
Stanic et al. (2009), who find a statistically significant 
difference in the £[oni] distributions of their optically- 
selected Revised Shapley-Ames sample of Seyfert Is and 
Seyfert 2s. 

4. EMISSION-LINE RATIO DIAGNOSTIC DIAGRAMS 

BPT emission-line ratio diagnostic diagrams can be 
used to infer whether the gas in a galaxy is primarily 
being heated by star formation or by radiation from an 
accretion disk around a central supermassive black hole. 
As discussed in Section 1, sources whose signal is dom- 
inated by accretion lie to the upper-right of those dom- 
inated by star formation. BPT diagrams are typically 
only applied to narrow-line active galaxies, because in 
BLAGNs the narrow lines are overwhelmed by emission 
from the broad-linc region. 

In Figure 3 we plot log( [OIII] /H/3) versus 
log([NII]/Ha) for the 21 sources in our 2 - 8 kcV 
sample for which we were able to determine the equiv- 
alent widths of all four emission lines (note that all 
21 have L x > 10 42 erg s" 1 ). Of the seven z < 0.5 
non-BLAGNs with Lx > 10 43 erg s _1 and measured 
[OIII] luminosities, only two appear in Figure 3. We note 
that this is a result of observational complications — sky 
lines interfering with our measurement of the H/3, Ha, 
or [Nil] lines or because the Ha/ [Nil] lines lie within 
the noise at the red edge of our spectral window. 
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The symbols grow in size and change in color with in- 
creasing X-ray luminosity. The dotted curve traces the 
Kewley et al. (2001) division between AGNs and extreme 
starbursts. It is a theoretical curve based on photoion- 
ization models for giant HII regions and a range of stellar 
population synthesis codes. The dashed curve traces the 
Kauffmann et al. (2003) division between AGN and star 
formers, which is based on an empirical separation of 
SDSS galaxies. 
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Fig. 3. — (a) A BPT diagram for the non-BLAGNs in our 2 - 
8 keV sample with z < 0.5. The symbols grow in size and change in 
color with increasing X-ray luminosity, as indicated in the legend. 
The dashed curve denotes the Kauffmann et al. (2003) division 
between AGNs (upper-right) and star forming galaxies (lower-left). 
The dotted curve denotes the Kewley et al. (2001) division between 
AGNs and extreme starbursts. 



Less than half (10/21) of the X-ray selected sources lie 
to the upper-right of the Kewley et al. (2001) curve and 
form a sequence similar to that of the SDSS Seyfert 2s, 
emerging from the HII region sequence and extending 
to the upper right-hand side of the diagram. Hereafter, 
we refer to these sources as 'BPT AGN-dominated non- 
BLAGNs'. We also refer to the sources that lie to the 
lower-left of the Kauffmann et al. (2003) division as 'BPT 
SF-dominated non-BLAGNs' and those that lie between 
the two curves as 'BPT transition non-BLAGNs'. We 
note that of the three BPT SF-dominated non-BLAGNs 
with 10 42 < L x < 5 x 10 42 erg s _1 (black triangles), 
none have L x < 2.5 x 10 42 erg s" 1 . A SFR of > 500 M 
is required to produce this level of X-ray luminosity (see 
Equation 2). Based on the space density of such strongly 
starbursting sources, we do not expect any in our z < 0.5 
2 — 8 keV sample. 

1 AGN lie in the transition 



Since two L x > 10 43 erg s 

vt2 



region and two 5 x 10 42 < L x < 10 43 erg s^ 1 AGNs lie 
below the Kauffmann et al. (2003) line, it is clear that 
the optical diagnostic diagram is not able to classify all of 
the X-ray selected non-BLAGNs with line ratio measure- 
ments as AGNs. Given that we are working from a high 
X-ray luminosity sample, this seems rather surprising. 
However, since it has been shown by many authors that 
there is considerable spread in £[oiirj/£jc (e.g., Alonso- 
Herrero et al. 1997; Bassani et al. 1999; Guainazzi et al. 
2005; Heckman et al. 2005; Ptak et al. 2006; Netzer et al. 
2006; Cocchia et al. 2007; Melendez et al. 2008; Diamond- 
Stanic et al. 2009; La Massa et al. 2009), it is worth ex- 



TABLE 1 
l°g(£[Oin]/ i x) 



Mean a 



BPT AGN-dom non-BLAGNs (z < 0.5) -1.06 (-0.36) 0.6 (0.6) 

Non-BLAGNs (z < 0.5) b -1.46 (-0.72) 0.6 (0.7) 

Non-BLAGNs (z < 0.85) b -1.66 (-0.97) 0.5 (0.6) 

BLAGNs (z < 0.5) -1.85 (-1.27) 0.5 (0.8) 

BLAGNs (z < 0.85) -1.76 (-1.14) 0.5 (0.7) 

a The values in parentheses are the £[oin],corr results. 
b Mean and a do not include absorber upper limit values. 

ploring whether this might affect how much overlap there 
is between the two AGN selections. 
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(b) Extinction Corrected L[0lll] 



QZiL_ 

f 1 1 1 1 1 1 in 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ifm 1 1 1 1 1 1 1 1 1 



-3 -2-10 1 -3 -2-10 1 
log (L[Olll]/Lx) log (L[Olll]corr/Lx) 

Fig. 4. — (a) Histogram of the log(L[oni]/^x ) for the sources in 
our 2 — 8 keV sample (solid line — z < 0.85 BLAGNs; dashed line — 
z < 0.5 BPT AGN-dominated non-BLAGNs). The arrow marks 
the mean value of the Heckman et al. (2005) extinction-uncorrected 
[OIII]-bright optically selected sample of Seyfert Is. (b) Same as 
(a) but for log(L[oill],corr/^x)- The arrow marks the mean value 
of the Mulchaey et al. (1994) extinction-corrected optically selected 
sample of Seyfert Is. 



5. THE RELATION BETWEEN L [0 iii] AND L x 
In Figure 4 we show the distribution of (a) 
log(£[om]/£x) and (b) log(L[oiirj, C orr/£x) for the BPT 
AGN-dominated non-BLAGNs (dashed line). They 
cover two orders in magnitude, confirming previous re- 
sults that showed a very large spread in this ratio. We 
also show the distribution for the z < 0.85 BLAGNs 
(solid line). 

In Table 1 we give the mean and standard deviation in 
\og(L [OUI] /L x ) and \og(L [om] cm /L x ) for the z < 0.5 
BPT AGN-dominated non-BLAGNs, for the z < 0.5 and 
z < 0.85 non-BLAGNs, and for the z < 0.5 and z < 0.85 
BLAGNs. In determining the mean values for the z < 0.5 
and z < 0.85 non-BLAGNs, we do not include the upper 
limit values for the absorbers (see Section 3.1). We find 
that the mean values for the z < 0.85 BLAGNs agree 
with those from the optically selected Heckman et al. 
(2005) extinction-uncorrected and Mulchaey et al. (1994) 
extinction-corrected Seyfert 1 samples (as indicated by 
the arrows in Figure 4). 

In Figure 5 we plot (a) log(L[oin]/Lx) versus log L x 
and (b) log(L[ iii], C orr/^ J '(:) versus log L x for the non- 
BLAGNs in our 2 — 8 keV sample. The black small tri- 
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(b) Extinction Corrected L[Olll] J 




42.0 42.5 43.0 43.5 44.0 

LOG (REST-FRAME 2-8 keV LUMINOSITY [erg s"']) 

Fig. 5.— (a) log(Z/[ in]/£x) vs. log L x for the L x > 
10 42 erg s -1 non-BLAGNs in our 2 — 8 keV sample (black small 
triangles — z < 0.5 non-BLAGNs lacking BPT line ratio measure- 
ments (excluding absorbers); black large triangles — z < 0.5 BPT 
AGN-dominated non-BLAGNs; stars — z < 0.5 BPT SF-dominated 
non-BLAGNs; circles — z < 0.5 BPT transition non-BLAGNs; blue 
small triangles — 0.5 < z < 085 non-BLAGNs). The upper limits 
indicate the absorbers in our 2 — 8 keV sample (see Section 3). The 
L[oiii] values are not extinction-corrected. The black (green) solid 
line shows the linear best fit to the z < 0.5 BPT AGN-dominated 
non-BLAGNs (the z < 0.5 BPT AGN-dominated sources plus the 
0.5 < z < 0.85 non-BLAGNs). The dashed curves show the 95% 
confidence intervals. We do not include the absorber upper limit 
values in our best fit. (b) Same as (a) but for £[oiii] extinction- 
corrected values. 

angles denote z < 0.5 non-BLAGNs lacking BPT line ra- 
tio measurements (excluding absorbers), the black large 
triangles denote z < 0.5 BPT AGN-dominated non- 
BLAGNs, the stars denote z < 0.5 BPT SF-dominated 
non-BLAGNs, and the blue small triangles denote 0.5 < 
Z < 0.85 non-BLAGNs. The upper limits indicate the 
absorbers in our 2 — 8 keV sample. 

Using a sample of 52 X-ray selected type 2 AGNs from 
the Chandra Deep Field-South and HELLAS2XMM, 
Netzer et al. (2006) found that log(L[oni]/-£<x) de- 
creases with increasing X-ray luminosity [linear best fit: 
log(£[oni]/£x)= -(0.42 ± 0.07) \ogL x + (16.5 ± 2.9)]. 
However, only eight of their sources have direct [OIII] 
measurements. For the others, they determined [OIII] 
using transformations from other available narrow emis- 
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LOG (Rest-frome 2-8 keV Luminosity [ergs s"']) 

Fig. 6.— log(I/[ OII T]/ix) vs. log L x for the L x > 10 42 erg s" 1 
BLAGNs with z < 0.5 (black squares) and 0.5 < z < 0.85 (blue 
squares) in our 2 — 8 keV sample. The £[om] values are not 
extinction-corrected. The black (green) solid line shows the lin- 
ear best fit to the z < 0.5 (z < 0.85) BLAGNs. The dashed curves 
show the 95% confidence intervals. 

sion lines. For example, for 32 of their sources they mea- 
sured the [Oil] emission line flux and then converted it 
to [OIII] using the mean [OIII] /[Oil] flux ratio from Za- 
kamska et al. (2003). They do not take into account 
any luminosity dependence in this ratio, nor is it known 
whether the Zakamska et al. (2003) mean value applies 
to an X-ray selected sample. 

In Figure 5 the black (green) solid lines show the best 
fit to the z < 0.5 BPT AGN-dominated non-BLAGNs 
(z < 0.5 BPT AGN-dominated non-BLAGNs plus the 
0.5 < z < 0.85 non-BLAGNs). The dashed curves pro- 
vide the 95% confidence intervals. We do not include the 
absorber upper limit values in our best fit. In Table 2 we 
give these best fit and confidence values. 

Using the Spearman Rank analysis we find a 1.6 a sig- 
nificance against a null result for the extinction-corrected 
z < 0.5 BPT AGN-dominated non-BLAGNs, which in- 
dicates little or no luminosity dependence. The Pearson 
coefficient of —0.6 indicates a weak anti-correlation (a 
Pearson coefficient of —1 indicates an anti-correlation, 
a coefficient of indicates a null result, and a coeffi- 
cient of 1 indicates a positive correlation). If we in- 
clude the 0.5 < z < 0.85 non-BLAGNs, then the sig- 
nificance against a null result goes up to 3.3 a, although 
the Pearson coefficient remains at —0.5. The slope for 
the z < 0.5 BPT AGN-dominated non-BLAGNs plus 
the 0.5 < z < 0.85 non-BLAGNs is consistent within the 
errors with the Netzer et al. (2006) results. 

In Figure 6 we plot log(L[oiir]/£x) versus log Lx for 
the BLAGNs with z < 0.5 (black squares) and 0.5 < 
z < 0.85 (blue squares) in our 2 — 8 keV sample. The 
black (green) solid line shows the best fit to the z < 0.5 
(z < 0.85) BLAGNs. The dashed curves provide the 95% 
confidence intervals. In Table 2 we give these best fit 
and confidence values. The slopes are consistent within 
the errors with the non-BLAGNs. Using the Spearman 
Rank analysis we find a 2.1 a significance against a null 
result for our extinction-corrected z < 0.85 BLAGNs and 
a Pearson coefficient of -0.3, indicative of little or no 
luminosity dependence. 
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TABLE 2 

Best fits to the luminosity dependent ratio of I/roini/i'X 



^g(L[oni]/Lx) = m \og(L x ) + b 
m b 
BPT AGN-dom non-BLAGNs (z < 0.5): 
£[0111] 

^[OrnLcorr 

BPT AGN-dom plus 0.5 < z < 0.85 non-BLAGNs 



-1.77 ±0.6 
-1.48 ±0.7 



74.1 ±26.5 
62.7 ±31.9 



S-cr a 
2.1 

1.6 



P-cocff b 



^[OIII] 



-0.54 ±0.1 
-0.60 ±0.1 



21.9 ±5.7 
23.0 ±5.9 



3.3 
3.3 



-0.7 
-0.6 

-0.5 
-0.5 



BLAGNs (z < 0.5): 
£[Olll] -0-45 ± 0.4 

^[OIII],corr -0.55 ±0.4 

BLAGNs (z < 0.85): 



17.5 ±16.6 
22.4 ±16.7 



0.1 

1.3 



-0.3 
-0.3 



J [OIII] 
^[OIII], 



-0.32 ±0.2 
-0.40 ±0.2 



12.1 ± 9.9 
16.1 ± 9.8 



1.5 
2.1 



-0.2 
-0.3 



a Spearman-Rank analysis. 
b Pearson coefficient. 

c Best fit does not include absorber upper limit values. 



6. DISCUSSION 

From our sample, which is a high Lx sample, we 
saw that there is substantial spread in L[ iii]/-^x at a 
given Lx for the z < 0.5 BPT AGN-dominated non- 
BLAGNs (see black large triangles in Figure 5), which 
means the sources can scatter all the way from the AGN 
portion of the BPT diagram to the star formation por- 
tion. This means that AGNs at the low end of the sample 
in L[oiii]/^x could be misidentified as star formers with 
the BPT diagram. 

In Figure 7 we replot the BPT diagram of Figure 3, this 
time using the color and size of the symbols to indicate an 
increasing log(L[oiii]/£x) rather than an increasing Lx- 
It does appear that there is a trend in log(L[ ni]/ix), 
with the z < 0.5 BPT SF-dominated non-BLAGNs ex- 
hibiting the lowest values (—3 <log(L[oin]/£x)< —2). 
This is consistent with the idea that the [OIII] emission 
from some AGNs is low. This dispersion in the amount 
of [OIII] light produced by an AGN also implies that we 
will not be able to tell the strength of the AGN relative 
to the star formation as we move along the wing, as some 
modelers have proposed (e.g., Kewley et al. 2006). In- 
stead the location along the wing will be a mixture of 
the AGN strength and how much L[om] there is for the 
AGN strength. 

One possible explanation for the low [OIII] emission 
in some AGNs is extinction from the host galaxy. In- 
deed, in Figure 5 we see that one of the five BPT SF- 
dominated non-BLAGNs (stars) significantly increases 
its log(L[oni]/ix) value as a result of the extinction cor- 
rection. However, the other four retain significantly lower 
log(L[ iii]/ix) values than the mean for the BPT AGN- 
dominated non-BLAGNs. We note that all five BPT 
SF-dominated non-BLAGNs have reliable Ha and H/3 
flux measurements and therefore unique extinction cor- 
rections. The small extinction corrections suggest that 
host galaxy obscuration is not the dominant reason for 
the low [OIII] emission in the majority of our BPT SF- 
dominated non-BLAGNs. 

Figure 8 also provides support for the idea that apply- 
ing extinction corrections does not help in reducing the 
scatter in the ratio between L[om] and Lx- In this fig- 
ure we show the distribution of (a) log(L[oiii]/^x) and 
(b) log(L[oni],corr/ix) for the z < 0.5 BLAGNs in our 




-1.0 -0.5 0.0 

LOG ([Nll]/Ho) 



0.5 



Fig. 7.— A BPT diagram for the non-BLAGNs in our 2 - 8 keV 
sample with z < 0.5. The symbols grow in size and change in color 
with increasing i[oin]/^Xi as indicated in the legend. The dashed 
curve denotes the Kauffmann et al. (2003) division between AGNs 
(upper-right) and star forming galaxies (lower- left). The dotted 
curve denotes the Kewley et al. (2001) division between AGNs and 
extreme starbursts. 
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Fig. 8. — (a) Histogram of the log(L[oill]/^x) for the z < 0.5 
sources in our 2 — 8 keV sample with reliable He* and H/3 flux mea- 
surements (solid line — BLAGNs; dashed line — non-BLAGNs). 
The arrow marks the mean value of the Heckman et al. (2005) 
extinction-uncorrected [OIII]-bright optically selected sample of 
Seyfert Is. (b) Same as (a) but for log(L[oill], corr/^x)- The ar- 
row marks the mean value of the Mulchaey et al. (1994) extinction- 
corrected optically selected sample of Seyfert Is. 

2 — 8 keV sample (solid line). We also show the distri- 
bution for our z < 0.5 non-BLAGNs (dashed line). We 
only include sources with reliable Ha and H/3 flux mea- 
surements. While the numbers are small, we see that 
correcting for extinction does not reduce the dispersion. 

It is possible that the extinction corrections determined 
from the Balmer line ratios in Section 3.2 are not totally 
correct for the [OIII] emission. For example, there may 
be additional contributions to the Balmer emission lines 
from star formation further out in the galaxy, so if the 
dust were not mixed uniformly throughout, the extinc- 
tion corrections measured from the Balmer line ratios 
might be less than what they should be for the narrow- 
line region. However, at least the extinction corrections 
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from the Balmer line ratios give some impression of how 
extinction might alter things, and we see that it does 
not solve the problem. Another issue is whether X-ray 
absorption due to an obscuring torus around the AGN 
could be affecting the luminosity ratio and whether cor- 
recting Lx for absorption by the torus would reduce the 
dispersion. 

We address these issues in Figure 9, where we show 
log(L[ iii]/ix) versus log L x for the BLAGNs in our 2 — 
8 keV sample as well as for the BLAGNs in other X-ray 
selected surveys from the literature. We see that a similar 
high dispersion is also observed for the BLAGNs, even 
though X-ray absorption would be small and extinction 
from the galaxy should not be very severe if the AGN is 
aligned with the disk of the galaxy. 
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LOG (Rest-frome 2-8 keV Luminosity [ergs s"']) 

Fig. 9. — log(L[om]/ix) vs. log L x for the L x > 10 42 erg s" 1 
BLAGNs with z < 0.85 in our 2 — 8 keV sample (black squares). 
The horizontal dotted line denotes the mean (log(L[oin]/-^x )) = 
— 1.76 for these sources. The £[om] values are not extinction- 
corrected. The green squares show the Heckman et al. (2005) sam- 
ple of z < 0.2 Seyfert Is from the RXTE slew survey (Sazonov & 
Revnivtsev 2004). The red squares show the Cocchia et al. (2007) 
Seyfert 1 sample from the HELLAS2XMM survey. 

We conclude that sometimes [OIII] will diagnose the 
presence of an AGN and sometimes it will not, and ex- 
tinction is not the primary factor. Furthermore, we find 
that the dispersion in L\p\\{\/ L x is larger than the dis- 
persion in L pt/Lx, which suggests that £[oin] is par- 
ticularly problematic. We characterize the optical con- 
tinuum luminosity, L opt , in our sources by their rest- 
frame g-band luminosity, L g . We determine L g by trans- 
forming the observed g,r,i, z, J, H, K s OPTX photom- 
etry (Trouille et al. 2008) into g-band flux at z = 
using kcorrect v4_l_4 (Blanton & Roweis 2007). Since 
the contribution from the emission lines in AGNs is so 
low compared to the continuum luminosity (see Heck- 
man et al. 2004 in which (L opt /L [om] ) ~ 320 for SDSS 
type 1 AGNs), L g provides a reasonable proxy for the 
optical continuum luminosity in these sources. We find 
a mean log(i op t/ 'Lx) for our BLAGNs of 0.94, with a 
dispersion of 0.36. This is smaller than the dispersion in 
log(L [0 iii]/ix) (<x = 0.5 for our BLAGNs, see Table 1). 

Instead we suggest that there is substantial complexity 
in the structure of the narrow-line region, which causes 
many ionizing photons from the AGN not to be absorbed. 



This complexity could include a low covering fraction of 
material, the configuration of the HII regions, the pres- 
ence of density-bounded HII regions, and wind venting. 
Hence, L\Qxa/L x would be low, and one would perceive 
the source as a star former rather than as an AGN. 

Several groups have compared the hard X-ray and 
[OIII] luminosity functions (LFs) for AGNs in the local 
universe (e.g., Heckman et al. 2005; Georgantopoulos & 
Akylas 2010) to try to determine which selection of AGN 
is most efficient. As we have shown, the large dispersion 
in L[ iii]/ix means that this transformation cannot be 
done to better than a factor of ~ 3 (one sigma) for an 
individual object. However, even though we cannot be 
confident about the tranformation, it is still interesting 
to explore whether the two LFs have, on average, the 
same shape. 

In Figure 10(a) we show the SDSS [OIII] LF for z < 
0.15 Seyfert 2s (Hao et al. 2005, black curve) and z < 0.3 
type 2 quasars (Reyes et al. 2008, blue curve). Following 
Reyes et al. (2008), we shift the Hao et al. (2005) lumi- 
nosities by 0.14 dex to account for the difference between 
the spectrophotometric calibration scales. Reyes et al. 

(2008) stress that their type 2 quasar LF is a lower limit 
due to their selection criteria, the details of the spec- 
troscopic target selection, and other factors. On this 
we show the Sazonov & Revnivtsev (2004) 3 — 20 keV 
LF based on their fit to the local RXTE sample (green 
curve). We also show the z = 0.1 (red solid curve; this 
is the volume-weighted redshift that corresponds to the 
redshift range of the SDSS Seyfert 2s) and z = 0.2 (red 
dashed curve; this is the volume-weighted redshift that 
corresponds to the redshift range of the SDSS type 2 
quasars) extrapolations of the Yencho et al. (2009) best 
fit to their distant Chandra plus local SWIFT BAT 
2 — 8 keV sample over the redshift interval < z < 1.2 for 
an independent luminosity and density evolution (ILDE) 
model. We have transformed the X-ray LFs to the [OIII] 
LFs using (log(i[oni]/ix)) = —1-06 (from the z < 0.5 
BPT AGN-dominated non-BLAGNs). The gray shading 
shows the 3 a uncertainties on the transformation. Since 
neither Hao et al. (2005) nor Reyes et al. (2008) apply a 
reddening correction to their [OIII] luminosities, we use 
our extinction-uncorrected transformation. 

The bottom panel shows the transformed Yencho et al. 

(2009) LFs divided by the SDSS LFs. For L [om] < 
10 42 erg s"\ we divide the Yencho et al. (2009) z = 0.1 
LF by the Hao et al. (2005) z < 0.15 LF. For i [OIII] > 
10 42 erg s~\ we divide the Yencho et al. (2009) z = 0.2 
LF by the Reyes et al. (2008) z < 0.3 LF. 

In Figures 10(b) and (c) we instead did the trans- 
formations using (log(L[om]/^x)) = —1-46 (from the 
z < 0.5 non-BLAGNs, excluding the upper limit values 
for the absorbers) and (log(Zroin]/£x)) = —1-85 (from 
the z < 0.5 BLAGNs). 

Despite the fact that there are strong selection effects 
in both the [OIII] and hard X-ray LFs (particularly in 
the way the [OIII] LFs were constructed, since they con- 
tain only Seyfert 2s, and due to the X-ray selection miss- 
ing heavily X-ray absorbed, Compton-thick sources), we 
see reasonable agreement. Thus, unlike Heckman et al. 
(2005), we do not see any evidence for a substantial in- 
crease in the [OIII] LF relative to the X-ray LF (see also 
Georgantopoulos & Akylas 2010). 
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Fig. 10.— (a) Comparison of the X-ray and [OIII] LFs (black curve^Hao et al. (2005) [OIII] LF derived from the z < 0.15 SDSS 
Scyfert 2s; blue curve — Reyes et al. (2008) [OIII] LF derived from the z < 0.3 SDSS type 2 quasars; green curve — Sazonov & Revnivtsev 
(2004) 3 - 20 keV local RXTE LF; red curve— Yencho et al. (2009) z = 0.1 extrapolation to their ILDE model best fit to the < z < 1.2 
rest-frame 2 — 8 keV distant -l-local LF; red dashed curve^ — same as the red curve but for a, z = 0.2 extrapolation). We aligned the X-ray LFs 
to the [OIII] LFs using the BPT AGN-dominated non-BLAGNs (log(L[oill]/^ / x)) = —1-06. The gray shading shows the 3 cr uncertainties 
for the transformed Yencho et al. (2009) z = 0.1 LF. Bottom panel shows the transformed Yencho et al. (2009) X-ray LFs divided by the 
SDSS LFs. (b) Same as (a) but using the z < 0.5 non-BLAGNs (log(L[oiii]/i'x)) = —1.46, which excludes the upper limit values for the 
absorbers, (c) Same as (a) but using the z < 0.5 BLAGNs (log(L[oill]/^'x)) = —1.85. 
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7. SUMMARY 

We have used our highly optically spectroscopically 
complete OPTX sample of 2-8 kcV X-ray selected AGNs 
to see how well an empirical optical emission-line diag- 
nostic diagram docs in classifying the X-ray sources as 
AGNs. We find that roughly 20% of the X-ray AGNs that 
can be put on these diagrams are not classified as AGNs 
by the optical diagnostic diagram, if we use the Kauff- 
mann et al. (2003) empirical division between AGN and 
star formers. If instead we use the Kewley et al. (2001) 
theoretical division, over half of our X-ray AGNs are not 
classified as AGNs. We show that both the BLAGNs and 
the non-BLAGNs in the OPTX sample exhibit a large 
(two orders in magnitude) dispersion in the ratio of the 
[OIII] to hard X-ray luminosities, even after the applica- 
tion of reddening corrections to the [OIII] luminosities. 
This confirms previous results from the literature. The 
large observed spread in the luminosity ratio leads us 
to postulate that the X-ray AGNs that are misidenti- 
fied by the optical diagnostic diagram as star formers 
have low [OIII] luminosities due to the complexity of the 
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